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q ■ Abstract 

O | We argue that recent NA49 results on multiparticle distributions 

■ and fluctuations, as a function of the number of participant nucleons, 

suggest that percolation plays an important role in particle production 
at high densities. 



D j Recentely, the NA49 collaboration has presented results, from the CERN/SPS 

at 158 A GeV, on multiplicity fluctuations or, to be more precise, on V(n) / (n) , 



X : V(n)Kn) S ( " 2) - ( " > '' 



V(n)l(n) S ^J^, (1) 

as a function of the number N part , of participant nucleons, from pp to PbPb 
collisions [1]. 

These data are very interesting for several reasons: 

1) They show evidence for universal behaviour: the experimental points 
in the plot V(n)/ (n) versus N part . fall into a unique curve (see Fig.l). 

2) The physics in the small N part _ limit (pp, N part _ — ► 2) and in the large 
N par t. limit (PbPb, N part , — > 2ApbPb) seems to be quite the same, as in both 
cases the quantity (1) approaches 1. The fluctuations are larger in the inter- 
mediate N part , region (see Fig.l). 
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Figure 1: Variance over average multiplicity, for negative particle production, 
as a function of the number of participants. The curves are from (16) (full 
curve, for NA49 data, dashed curve for RHIC, y/s = 200GeV). 

3) The (negative) particle distribution, in the low density and in the high 
density limits, is in fact a Poisson distribution (see Fig. 2), the distribution 
being wider than Poisson in the intermediate N part . region. 

In the framework of the string model with percolation [2], these results are 
natural. On one hand, percolation is a universal geometrical phenomenon, 
the properties depending essentially on the space dimension (dimension 2, 
impact parameter plane, in our case), and being controlled by the transverse 
density variable i], 



where r is the transverse radius of the string (r ~ 0.2/m), R the radius of 
the interaction area, and N s the average number of strings. The quantity 
(R/r) 2 is nothing but the interaction area in units of the string transverse 
area. As R and N s are functions of the number N part . of participants, N part , 
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Figure 2: Multiplicity Distributions, P( function of n_. The curves 

are Poisson (dashed lines) and Negative Binomial (full line). 
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similarly to r], becomes, at a given energy,a universal variable. 

On the other hand, in percolation [3], what matters is the fluctuation in 
the size of the clusters of strings: one starts, at low density (small N part _), 
from a situation where strings are isolated, at intermediate density one finds 
clusters of different sizes, and one ends up, at high density, above the perco- 
lation threshold, with a single large cluster. In both, low and high, density 
limits, for a fixed number N s of strings, fluctuations in cluster size vanish 
(see Fig. 3). In the simplest string model the particle distribution is Poisson 
(as observed in e + e~ and pp at low energy) and V/ (n) — > 1 in both, low and 
high, density limits (as seen in Figs.l and 2). 

In hadron-hadron and nucleus-nucleus collisions, during the collision strings 
are produced along the collision axis, and these strings may overlap and form 
clusters of different sizes. In the spirit of percolation theory, we shall assume 
that fluctuations in the number N of strings per cluster dominate over all 
the other fluctuations. 

Following [4] we write for the multiplicity distributions at a given value 
of r] (or N part ) 

P{n) = J W(X)p(n,Xn)dx , (3) 

where the integral is over the cluster distribution W(X), and p(n,Xn) is a 
convolution of Poisson distributions, n being the single string average mul- 
tiplicity. If W(X) is a gamma function, P(n) becomes a negative binomial 
(NB) distribution. From (3) we obtain 

(n) = (X)n , (4) 
1/fc _ (n 2 )-(n) 2 1 _ (X 2 ) ~ {X? 

1 W V) - {x? ' (5) 

k being the NB parameter. As, see [5] , 

X = N C N , (6) 

where N c is the average number of clusters and iV the number of strings in 
a given cluster, and the sum-rule, 

N s = N C (N) (7) 

is valid, we have, instead of (4) and (5), 
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Figure 3: Impact parameter percolation. For small densities (77 <C 1) and 
for large densities (77 ^> 1) there are no strong fluctuations in the number N 
of strings per cluster. For intermediate densities (77 ~ 1) iV-fluctuations are 
large. 

(n) = iV s n , (8) 

and 

As the NA49 experiment is at relatively low energy (s/s ~ 20GeV) we 
shall have contributions from valence strings (V) and from sea strings (5). 
We thus write, 

N v N s 
n = ^-f + jf-n s , (10) 

with 

N s = NY + N? , (11) 

with, as valence strings are longer, n v > n s . We then have, from (8) and 
(10), 

(n) = N^n v + N^n s , (12) 

and, from (1) and (9), 

(n) k 
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Figure 4: The parameter A; as a function of r/: the "coins in boxes" model 
(dashed curve) and the percolation model (full curve). 



^From the percolation arguments given above, (iV 2 ) — (iV) 2 — > as r\ — > 
and i] — > oo, which implies k(r)) — > oo (Poisson distribution) in the same 
limits. As (n), similarly to N s , is a monotomically increasing function of i] 
(and of energy and of N part .) it is clear that from (13) one expects a maximum 
of V(n)/{n) at same value of r\ and V(n)/(n) — > 1 at low and density as 
seen in data. 

An attempt was made in [5] to study the k dependence on 77 in pp{pp) 
collisions. Experimentally at low 77 (or energy) k is large and decreases with 
increasing 77. This is true in percolation and is true as well in the "coins-in- 
boxes" model utilized in [5]. For the fluctuation of the number N of coins 
per box one obtains in that model 

HV) = -1 n \ 1 V ( 14 ) 

v ' 1 - (1 + rj)exp(-rj) 

and the curve (14) is shown in Fig. (4) - dashed line. At low density 77 
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the "coins in boxes" model is similar to percolations: strings are isolated, 
coins sit in different boxes. At high density, as the boxes do not loose their 
identity, the equivalent to the appearence of a single big cluster does not 
occur: this means that fluctuations do not vanish at rj — > oo, as it happens 
with percolation. 

In order to have a model satisfying the fix N s percolation conditions 



(N 2 ) - (N)' 







T] — > 0, Tj — > OO 

we modified the "coins in boxes" model, (14), to write 



k{rj) = b 



M _ i 



1 - (l + br])e 



(15) 



(16) 



such that (16) agrees with (14) at low 77, but k{r\) in (16) increases much 
faster at high density (see full curve in Fig. (4). 

Before making a comparison with NA49 data, there are two questions to 
be addressed: 

i) F(rj) factor due to random colour summation 

When strings fuse in a cluster the effective colour charge is not just the 
sum of the colour charges of the individual strings [6]. In pratice, the effective 
number N of strings is reduced, [7], 



N — ► F(i])N , 



where 



F( V ) 



rj 



such that, instead of (12), 

(n) = F(rj) \N^n v + N? n s 



(17) 
(18) 

(12') 



but (13) remains unchanged. 

ii) The relation between 77 and N part , 

In the definition of 77, (2), what appears is not N parL but rather the average 
number N s of strings and the radius R of interaction. Making use of simple 
nuclear physics and multiple scattering arguments, one has [8] 



R ~ R p n\ /3 



(19) 
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N s ~ N p s N A f , (20) 

and 

~ 2iV A . (21) 

where R p is the nucleoli radius (~ 1/m), JV| is the (energy dependent) 
number of strings in pp collisions, at the same energy, and Na is given by 

N A = . (22) 

^From (2), (19), (20) and (22) we obtain for the relation between rj and 

Npart. i 

We turn now to the comparison with experiment: 
1- Multiplicity (n) 

Instead of using directly (12') we shall write (n)N A as a function of (n) p , 
making use of (12'), (19), (20), (21), (22) and (23), 

<»>- ^ ^"^-^f ■ < 24 » 



2 (n y - n s ) 



with 

\n v - n~ 

being an adjustable parameter, c > 0, decreasing with energy. 
In Fig. (5) we show (17) in comparison with NA49 data. 

2- Fluctuations: V(n)/(n) 

In Fig. (1) we show the comparison of the percolation model (16) with 
the NA49 data, and the prediction at y/s — 200GeV. 

The parameters used in both fits, Figs. (5) and (1) were: b = 1.68, 
c = 2.3, (n) p = 0.52 (from data), Nf = 3.5, (R p /r) ~ 5. 

In Figs. (5) and (1) the dashed lines represent our prediction (without 
phase space adjustments) for RHIC (y/s = 200GeV), with Nf = 7.5 and 
(n) p = 0.70 (see first paper in [2]). 
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Figure 5: a) The multiplicity as a function of N part _ in comparison with NA49 
data (full line). Prediction at RHIC, y/s = 200GeV (dashed curve); b) The 
ratio {ti)/Na as function of N part . (full line for NA49 data and dashed line 
for v/i = 200GeV). 
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In conclusion, we find that the recent NA49 results, regarding the mul- 
tiplicity distribution dependence on the number of participating nucleons 
are quite consistent with the impact parameter, percolation description of 
hadron-hadron and nucleus-nucleus collisions at high density and high en- 
ergy, as proposed in [4] and [5]. For related work see [9] and [10]. 
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10 



References: 



1. M. Gazdzicki et al, J. Phys. G30, 5701 (2004) (nucl-ex, 0403023). 

2. N. Armesto, M.A. Braun, E.G. Ferreiro and C. Pajares, Phys. Rev. 
Lett. 77, 3736 (1996); M. Nardi and H. Satz, Phys. Lett. B442, 14 
(1998); A. Rodrigues, R. Ugoccioni and J. Dias de Deus, Phys. Lett. 
B458 (1999) 402. 

3. D. Stauffer, Phys. Rep. 54, 2 (1979); D. Stauffer and A. Aharony, 
Introduction to Percolation theory, Taylor and Francis (1992). 

4. J. Dias de Deus, E.G. Ferreiro, C. Pajares and R. Ugoccioni, Euro 
Phys. Journal, to appear, (hep-ph/0304068); C. Pajares, Acta Phys. 
Pol. B35, 153 (2004). 

5. J. Dias de Deus, E.G. Ferreiro, C. Pajares and R. Ugoccioni, Phys. 
Lett. B601 (2004) 125. 

6. T.S. Biro, H.B. Nielsen and J. Knoll, Nucl. Phys. B245, 449 (1984). 

7. M.A. Braun, F. Del Moral and C. Pajares, Phys. Rev. C65, 02490 
(2002). 

8. J. Dias de Deus and R. Ugoccioni, Phys. Lett. B491, 253 (2000), Phys. 
Lett. B494, 53 (2000). 

9. P. Brogueira and J. Dias de Deus, Acta Phys. Polon. B36 (2005) 307. 

10. L. Cunqueiro, E.G. Ferreiro, F. del Moral and C. Pajares, |hep-ph/ 0505197 
(2005). 



11 



